Inflation and dark energy from three-forms 
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Three-forms can give rise to viable cosmological scenarios of inflation and dark energy with po- 
tentially observable signatures distinct from standard single scalar field models. In this study, the 
background dynamics and linear perturbations of self-interacting three-form cosmology are investi- 
gated. The phase space of cosmological solutions possesses (super)-inflating attractors and saddle 
points which can describe three-form driven inflation or dark energy. The quantum generation 
and the classical evolution of perturbations is considered. The scalar and tensor spectra from a 
three-form inflation and the impact from the presence of a three-form on matter perturbations 
are computed. Stability properties and equivalence of the model with alternative formulations are 
discussed. 
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I. INTRODUCTION 

Inflation is a successful explanation of many cosmolog- 
ical puzzles, and the current acceleration of the universe 
is a cosmological puzzle which yet lacks an explanation. 
Since Nordstrom |l| , scalar fields have been present in ex- 
tra dimensional and fundamental theories, and it is natu- 
ral to employ them to describe the energysources needed 
to generate inflation and dark energy 0, d, @, 0] j for 
recent reviews, see 18|,19|]. However, it is crucial to under- 
stand how strict are the theoretical and phenomenologi- 
cal limits on the role of higher spin fields in cosmology. 

Vector inflation has been considered recently, us- 
ing either time-like [Ul or space-like components. 
However, to naturally inflate, the vector needs a non- 
minimal coupling and seems to feature instabilities p^ . 
see however Il5l |. Effects on CMB, alternative sce- 
narios [ig'I and the perturbation generation [13, IH, fioj 
have been studied in these models. Vector field dark en- 
ergy might alleviate the coincidence 



problem 25, 2u, '27[ and introduce new effects on per- 
turbations [28, 29]. It can be shown that two- form in- 
flation resembles much the vector inflation, having the 
same possibilities and problems [s^, [3l|, [111- Spinor 
[33I [34 1 and Yang-Mills [H, [s^ flelds have been also 
explored. Kalb-Ramond forms with dilaton couplings 
has been considered in the frameworks of strin g co smol- 
ogy dynamics [13, [111 , pre-big bang cosmolo gy l39ll. uni- 
fied models of dark matter and dark energy [40, [41| and 
bouncing cosmology |42| . T wo- forms appear also in the 
asymmetric gravity [43l |44| as the antisymmetric contri- 
bution to the metric, and have been considered in cos- 
mology (45l . l46j . Recently a Chern-Simons type gravity 
was developed promoting the Levi-Civita symbol into a 
dynamical field [131 ■ In the two- measures field theory. 



the new measure of integration can be built from either 
four scalar fields or an independent dynamical three-form 
481. This theory has several cosmological implications 
49|. In a scale- invariant realization of the theory, a 
scalar potential acquires two fiat regions in such a way 
that both inflation and dark energy may emerge [HO, l5ll |. 
Forms, being intrinsically anisotropic, could also be rele- 
vant for a dynamical origin of the four large dimensions 
modelin g v iolation of the Lorenz invariance 
Js^JsqI. l60l , the observed CMB anomalies 




64, 651 or testable late-time anisotropic phe- 
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nomenology [66., £7,, |68|, l69|, |70| . 

In the present paper, our aim is to study the possible 
cosmological significance of three-forms. It was noticed 
in [rij that the four-form constructed from a three-form 
gauge potential generates a cosmological constant. Since 
then, this fact has been employed in discussions attempt- 
ing to explain the tiny (or vanishing) value of the cosmo- 
logical constant [zl, U^. Recently, we have proposed 
to consider the case of self-interacting gauge potential 
[3^ [tI ] . This breaks the gauge- invariance but the field 
becomes then dynamical. Then a single field inflation 
with an exit to radiation dominated era can be naturally 
generated, or alternatively, the three-form can act as pos- 
sibly transient dark energy at a late stage of the history 
of the universe (three-from induced potentials were dis- 
cussed in [75, 76]). Form fields appearing in string theory 
generically couple to branes and this way a potential term 
might be obtained. In the present study we conflne our 
investigations to the simple model with only a canonical 
field minimally coupled to Einstein gravity. A three-form 
generalization of vector (and scalar) infiation was intro- 
duced recently [s^ , based on an action involving nonmin- 
imal couplings in such a way that the equation of motion 
of the comoving field has exactly the Klein-Gordon form 
in FLRW spacetime. The study of gravitational waves 
in such model reveals an instability occurring at large 
values of the field 3l| while the spectrum of scalar per- 
turbations in small field inflation could be slightly red 
tilted and thus compatible with observations [7^ . In the 
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minimally coupled model that we consider here the equa- 
tion of motion can be also written in the Klein-Gordon 
form, but given an effective potential. 

The dual of the three-form is a scalar field. In the 
case of a non-quadratic potential, the kinetic term of the 
scalar field is noncanonical. Such a model then becomes 
equivalent to k-inflation, and has been analyzed before 
[tsI [roj . A non-quadratic dependence on the three- form 
Faraday term results in a self-coupling of the scalar field. 
For any non-minimal coupling, in particular nonminimal 
gravity couplings of the three-form, the duality with a 
scalar field breaks down. In fact, as we will show ex- 
plicitly, even for some fairly simple self-interactions of 
the three-form, no dual description can be established in 
terms of a scalar field. In any case, typically the scalar 
field description, if it exists, is opaque and intractable 
even if the model in three-form language is simple and 
intuitive. In fact most of the models admit a reformula- 
tion as vector models, and in some cases a presentation 
as dynamical four-form models is also possible. Thus the 
simple starting point we have opens new perspectives on 
several classes of form cosmologies. 

We will write down the basic equations in section |TT] 
and review some results of Ref. [74] to give an intuitive 
picture of the possible background dynamics. A conve- 
nient variable to describe is the comoving field X. An 
analogy between a scalar field and the comoving field X 
can be utilized to illustrate the behavior of the field, how- 
ever, since the kinetic term of the actual field X/a^ gov- 
erns some aspects of the dynamics, an alternative view- 
point is also necessary to fully understand the dynamics. 
In section IIIII we give a detailed account of the back- 
ground expansion dynamics. The phase space analysis of 
the system reveals three classes of fixed points, one corre- 
sponding to matter domination and two corresponding to 
the domination of the form field. The nature and stabil- 
ity of the latter two points, which are relevant for dark 
energy and inflation solutions, depends on the form of 
the potential. We consider exponentials, power-law and 
Ginzburg-Landau type potentials. 

In section ITVl we consider perturbations of three- form 
cosmology. Guided by the duality with a vector field, we 
parameterize the four degrees of freedom in the fluctua- 
tions of the three-form, two of them transforming as vec- 
tors under spatial rotations and two as scalars. It turns 
out that the vector-type fluctuations can be neglected 
about a FLRW background, whereas the scalar pertur- 
bations introduce several possible new effects in cosmol- 
ogy. We analyze the quantum generation of three-form 
perturbations during inflation using the standard tech- 
niques, and give the detailed general form of predictions 
for the amplitude of scalar and tensor fluctuation spec- 
trum and their spectral indices. The possible influence on 
matter inhomogeneities due to the presence of classical 
perturbations in a three-form in a matter-dominated uni- 
verse is also considered. It is found that depending on the 
sound speed of the X-component (which in turn depends 
on the chosen form of the potential) , there can be a range 



of scales where the linear growth of matter density per- 
turbation is affected by the three-form fluctuation. This 
effect can be encoded into an effective strength of grav- 
itational coupling of matter particles, which in general 
depends on time and length scale. 

Finally we discuss some formal aspects of the model. 
In section |Vl some manipulations of the action are per- 
formed in order to clarify the properties of the model. 
The degrees of freedom and their nature will be seen to 
depend strongly on the form of the self-interactions. For 
some specific cases, dualities and equivalences can be es- 
tablished with other form field models as stated above. 
The model is also discussed briefly in the wider context 
of general (quadratic) three-form actions. We conclude 
in section IVII stating a few central formulas we have ob- 
tained. Some details of the scalar field formulation in 
particular cases is given in the Appendix]^ and an alter- 
native viewpoint employing the dual vector is mentioned 
in appendix iBl 



II. EXPANSION DYNAMICS 

We shall focus on a canonical theory minimally coupled 
to Einstein gravity. The action for a three-form A can 
then be written 
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where — SttG. In this section, to avoid unnecessary 
and excessive use of indices [o^ , we introduce the follow- 
ing notations: squaring means contracting the indices in 
order, as — A^fijA'^^'^ , dotting means contracting the 
first index, as (V • ^)q/3 = V^A^ap, and circling means 
contracting all but the first index in order, [A o B)ij^i, = 
A^apB^,"^ and finally antisymmetrization is performed 
by square brackets, for example [Al\^„ = ^i^ti-f ~ 
Now the possible drawback seems to be that the valence 
of the objects we are dealing with is not explicitly seen; 
however, in this index-free notation most results auto- 
matically generalize for tensors of arbitrary valence. In 
this notation the F{A) is F = {n + 1)[VA], where A is 
a n-form, and thus F{A) generalizes the Faraday form 
appearing in Maxwell theory. The energy momentum 
tensor is 



r= ^F o F + &V{A^)Ao A~ g i^F^ + V{A^) 



(2) 



The action leads to the equations of motion 

\I -F = UV'{A^)A, (3) 

which implies, due to antisymmetry, the additional set of 
constraints 



V • V'{A^)A ^ 0. 



(4) 



We consider a flat FLRW cosmology described by the 
line element 



(5) 
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The nonzero components of most general three-form com- 
patible with this geometry are then given by 



(6) 



where we have considered, instead of the field A, the 
more convenient comoving field X, and i, j, k denote 
the spatial indices. The relation between the squared 
invariant and the comoving field X is then A^ = QX^ . 
We can thus consider the potential as a function of X^, 
as we will do in the following. 

The equation of motion of the field X is then: 



X = -3HX -Vx- iHX . 



(7) 



An overdot means derivative with respect to the cosmic 
time t. The background perfect fluid evolves with 



Pb = SjHpB , 



(8) 



where 7 = 1+pb/pb, and these equations are subject to 
the Friedmann constraint 



'L.(l(X + 3HXf + V{X) + pB 



(9) 



The other field equation follows also from this and the 
continuity equations as 



H = -—{VxX + jpb) 



(10) 



We can thus define energy density and pressure of the 
field as 



Px - ^iX + 3HX)^ + ViX)., 



(11) 



Px 



-{X + 'SHXf -V{X) + VxX . (12) 

(13) 



The equation of state parameter of the three-form, wx 
Px/px, can be written as 



Wx 



'1 + 



VxX 

PX 



(14) 



We thus see directly that whenever the potential or just 
its slope vanishes, the field is like a cosmological con- 
stant. Furthermore, whenever the slope of V{X) is neg- 
ative (positive) if X is positive (negative), the comoving 
field behaves as a phantom field. So the origin has some 
absolute meaning for this field, unlike in the case of a 
scalar. We also note that the equation of state is un- 
bounded from both up and below. 

Previously it was shown that one may predict the evo- 
lution of the system by considering the effective potential, 
defined by ^ 



offV 



V 



X 



3HX. 



(15) 
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FIG. 1: The potential (red, solid line) and the effective po- 
tential (blue, dashed line) for the potential V — {X^ — C^)^ , 

when C" = I < 1/1 and ps — 0. Units of «: = 1. 



C = L5 



10 



\ V={X^-C^f 








^ \ 1 
' » / 






>. ' > / 

\- ' / 
\ - - 1 7 

\ * / 
\ * / 
\ * / 

\ 7 

\ / 1 
^v—X 1 




Veff 


1 
1 
1 

1 • 
1: 

! 
1 
1 



-2 



-1 





X 



FIG. 2: The potential (red, solid line) and the effective po- 
tential (blue, dashed line) for the potential V — {X^ — C^)^ , 



when C = I > 



and Pb ~ 0. Units of «: = 1. 



We illustrate the form of this potential in Figs. [T] and 
[5] for two new cases when pB = 0. For these poten- 
tials, the positions of the minima depend on the precise 
value of C. When C < \/2/3, the local minima are 
at X = ±C and when C > •\/2/3 the minima are at 
C — ± a/2/3. This suggests that the field has different 
dynamics given a choice of C. In particul ar it seems that 
the late time value of X is ±C if C < ^2/3 and X ap- 
proaches if C > •\/2/3. We will see in the Section 
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mil that these first impressions are indeed correct. One 
notes that there are places where the slope of the effec- 
tive potential is downwards while the bare potential is 
increasing. In such situations the field, as it rolls down 
the effective potential, climbs up its bare potential. The 
analogy with the scalar field holds here in the sense that 
such cases indeed correspond to a phantom-like expan- 
sion, where the (effective) equation of state of the field is 
more negative than minus unity. 

We illustrate these behaviors also with numerical so- 
lutions. In Figs. ([3|) and ([4]), we show the evolution of 
the field, energy density and equation of state for the 
Landau-Ginzburg potential with C = 0.5 and C = 1.5, 
respectively. We use this potential because it illustrates 
most features of the dynamics. We will see that depend- 
ing on the magnitude of C, two qualitatively different 
cases emerge. Though the early evolutions are similar 
while the background is dominant and the potential is 
roughly V ~ X^, the late time evolutions diverge. This 
is to be expected because, as we have seen in Figs. (P) 
and the effective potentials are different. In particu- 
lar, we se e tha t the field settles in X = C m the first case 
as C < ■y/2/3 and the late time oscillations result in the 
equation of state parameter to cross wx = — 1 at each 
oscillation. This can be understood by noticing that the 
field is transiting between positive and negative values of 
the slope of the potential at the minimum. For the second 
case C > -^2/3, the field cannot reach its minimum as it 
is constrained to X < for positive velocity. In par- 

ticular, we see that X cannot be displaced further than 
X — ±-y/2/3 for X = 0, as this saturates the Friedmann 
equation ^ and larger field values could be reachable 
only in the presence of negative energy sources. Conse- 
quently the equation of state approaches asymptotically 
w = —1 from below in the second case where C is larger 
than the critical value. 

In Figs. [3] and [5] we see that at early times the evolu- 
tions are identical and initially are described by a track- 
ing behavior followed by a constant energy density of the 
field X and only at late times the evolutions diverge from 
one another. The earlier history of these evolutions cor- 
responds to power-law potential regime of the potential, 
which was discussed in detail in Ref. ^tH, but we make 
here brief comments from the dual point of view. 

It was shown that when the background fluid is dom- 
inant, and the one independent component of the three- 
form Aijk is constant, the field will be tracking with 
n = 7, where n is given by the shape of the potential, 
V{X) = FqX". This period is ended when m^y^ « 
V{X). It was found that the number of e- folds elapsed 
is then rather accurately given by 



N, 



(16) 



where now is the initial value of y = k{X' + 3X)/\/6 
and Vi/Hf is the initial ratio of the potential and the 
Hubble rate and prime means differentiation with respect 
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FIG. 3: Cosmological evolution as a function of the e-folding 
time A'^ = Ina for the Landau-Ginzburg potential when 
C = 0.5 (units of K = 1). In the upper left panel we see 
the field going through the unstable fixed point at X = -^2/3 
to oscillate around the stable minimum a.t X = C. The upper 
right panel shows the corresponding values of the potential. 
The lower left panel depicts the energy densities in a logarith- 
mic scale; the dashed (red) line is matter and the solid (blue) 
line is the three-form. A brief tracking period is included in 
the figure, followed by two stages: the unstable fixed point 
and the oscillations around the stable minimum. The lower 
right panel shows the equation of state of X, exhibiting strong 
oscillations during the settling to the minimum. 



to 



N = \na. 



(17) 



The numerical solutions shown here also follow this be- 
havior now corresponding to the case n = A. 

One may consider the tracking property in light of the 
dual description as a scalar field which we discuss in more 
generality in section |VBj Now such a dual description as 
a scalar field exists, since during the tracking phase the 
kinetic term is negligible and the potential can be approx- 
imated by a power-law. Then this dominating potential 
term represents the kinetic term of the scalar field. It can 
be easily shown that a k-essence Lagrangian C ^ {dcj))'^^ 
behaves as a perfect fiuid with the equation of state given 
by = l/(2p — 1). During the tracking phase described 
above, the three-form indeed is (at least approximately) 
equivalent to such a k-essence field. One may deduce 
from the action (|130[) that a power- law potential V ^ 
turns into the kinetic Lagrangian with 2p = n/(n ~ 1), 
thus assuming the scaling of energy density we have ob- 
tained. 

Canonical quintessence is also known to possess the 
tracking property in some cases. One might thus be 
curious if this fact could be used to construct tracking 
three-form models. It turns out to be the case, and to 
produce noncanonical three- forms models. Among the 
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FIG. 4: Cosmological evolution as a function of the e-folding 
time A'' = Ina for the Landau-Ginzburg potential when 
C = 1.5 (units of K = 1). In the upper left panel we see the 
field going through the unstable fi xed point at X = to settle 
to the stable minimum at X = ^2/3. The upper right panel 
shows the corresponding values of the potential. The lower 
left panel depicts the energy densities in a logarithmic scale; 
the dashed (red) line is matter and the solid (blue) line is the 
three-form. A brief tracking period is included in the figure, 
followed by two stages: the unstable fixed point and the ap- 
proach to the late time attractor X — ^2/3. The lower right 
panel shows the equation of state of X, exhibiting a phantom 
behavior during the approach to the stable minimum. 



simplest examples of a quintessence with a tracking at- 
tractor is an inverse power-law potential Vi^cj)) = Vo(t>P. 
It is known to approximately track the background den- 
sity rather independently of the initial conditions. The 
power-law form of the potential translates into a power- 
law kinetic term of a three- form 4", 



density can then scale as a power of the scale factor, 
given the initial condition that the (only independent) 
component of the field, Aijk is a constant, meaning that 
the comoving field scales as X ^ a~^. 

Before settling into the minimum, the field turns 
around the potential and start climbing it. The value 
of Nt when this happens is given by 



1 , (IB 

3(1 + 7/2) ""{lA 



(21) 



where A = ^/2/3yi/{l + '-f /2), and B = Xi-A, where Xi 
is the initial value of the field. To have scaling behaviour 
for many e-folds, one may consider tiny A at huge Xi. 

Let us then consider the transient acceleration sce- 
nario. There the fie ld s lows down never reaching the 
critical value X = \/2/3, the amount of inflation Ne is 
eventually given by the initial value of the field Xi at the 
beginning of inflation near the critical point and again 
depends on the slope shape. 



KXi 



± 



1 



9n 1 + 2Ne 



1/2 



(22) 



and clearly, the slow-roll condition on the velocity of the 
field (|96p must be satisfied. We observe that for larger 
values of A^e then \KXi\ must be closer to \KXi \ « a/2/3. 

We also showed in Ref. [t^ that the oscillations of the 
field, when it settles to the minimum, follow an averaged 
behavior which depends on the shape of the potential. 
For a power law potential the result is, surprisingly, the 
same as for a canonical scalar field 1801, 



n- 2 



(23) 



VopP 



(18) 



where the kinetic term is given by 



The exponential potential is known to be the special case 
possessing the scaling, or "exact tracking" property. In 
section lVBI we will describe how the three-form can some- 
times be written as quintessence; by going the other way 
around one may find that a quintessence model specified 
by minimal coupling and V{(p) = Voe"^*^ can be recast 
into the three-form model 



In- 



- J_^2 

18 



(20) 



In the present study we however confine to the canonic 
kinetic term. As shown above, the three-form energy 



Thus for n = 2 the fields behaves as dust, (wx) = 
and for n = 4 it mimics radiation, (wx) = 1/3 [9y|. We 
mention in passing that this result applies also for oscil- 
lating k-essence with mass term and a power-law kinetic 
term due to the duality mentioned above, so the model 
(19) £ ~ (a</>)2p-hm202 oscillates like {w^) = (p- l)/(3p- 1). 

This concludes our qualitative review of the possible 
sequences of cosmological epochs, and next we turn into 
more rigorous phase space analysis and observable pre- 
dictions. 



III. PHASE SPACE ANALYSIS 

In this section we will put on a more solid and formal 
ground the considerations of the previous section on the 
late time dynamics of the system and its stability. We 
start by rewriting the equations of motion in the form of 
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TABLE I: 
The critical points in the system. 

a system of first order differential equations 




(24) 



xy 



(7 + A(x)(l - ~ w^) X - 7w^) , (26) 



where we have defined 
X = kX , 



K ,,,, 2 



y=^(X' + 3X), 



.,2 _ 



, (27) 



(28) 



A(a;) is, therefore, a function of X. The quantity z was 
eliminated from the equations of motion by applying the 
Friedmann constraint 



2 I 2 

y +z 



w 



1. 



(29) 



The system has three critical points which 

are described in Table [J 

A: a; = 0, J/ = 0, w ±1, for any A. ft corresponds 
to the background dominated solution. At late time the 
ratio H / H'^ approaches —87/2. ff 7 is a constant, the 
eigenvalues are (—3,37/2,87), hence, it is an unstable 
critical point. 

B: a; = ±-\/2/3, y — ±f, w — Q, for any A. This is a 
critical point that does not exists for the standard scalar 
field and that results from the extra X dependent terms 
in the equation of motion and in the definition of the en- 
ergy density and pressure. When approaching this fixed 
point, H approaches a constant at late times, however, 
keeps increasing, therefore, the effective equation of 
state parameter of the field X approaches —1 from be- 
low. The eigenvalues are (—3, 0, —87/2) and because one 
of the eigenvalues is zero, we cannot infer anything about 
the nature of the critical point from the linear analysis. 
We need to consider specific potentials and go to nonlin- 
ear order. The eigenvector corresponding to the vanish- 
ing eigenvalue reads (^^2/3, 1, 0), therefore, when going 
to higher order we study the stability of perturbations 
along the zero eigenvalue direction 5r ~ ^J2/'i5x + 6y, 
for which we get 



for n > 1 and is the coefficient resulting from ex- 
panding equations Eqs. ([M)) and (P5|) to nth order and 
using 5x = \/6<5r/5 and 5y — 3(5r/5, such that /i^^^ — 1. 
The general solution to this equation is 

5r = 5ro [I- Sr"^-^ {n - I) ^("^ iV 1 



(31) 



In order for a negative initial perturbation [Sr^ < 0) to 
decay one must have /x*^"^ > if n is even and /i'"-* < if 
n is odd. For a positive perturbation it suffices to have 
regardless of the value of n. 

Since \y \ must be less than unity, there can only be neg- 
ative perturbations along the r direction about the fixed 
point {x,y,w) — (■\/2/8, 1,0) and positive perturbations 
about the fixed point {x,y,w) — (— •\/2/3, —1, 0), thus, 
for the perturbations to decay it is required that /i^"-* 
must be positive for fixed point {x,y,w) = (■\/2/3, 1, 0) 
if n is even and negative if n is odd and /j,*^"-' negative 
for fixed point {x.y^w) = (— •\/2/3, — 1, 0). 

C: X — Xoxtj y = ^/S/ 2 a^ext : w = where Xext corre- 
sponds to the value of x at the extrema of the potential, 
i.e. where A = 0. In this case, becomes constant 
and H vanishes at late times. The stability of these fixed 
points, is therefore, strongly dependent on the specific 
form of the potential. 

We shall now look at particular examples to illustrate 
the significance and stability of the fixed points just de- 
scribed. 



1. V = exp{-l3X) 

Because the potential, being a function of the invariant 
A^, must depend explicitly on X^ instead of AT itself, we 
are dealing with symmetric potentials. This potential is 
not X^ dependent and therefore should only be seen as 
an example to be compared with the standard scalar field 
dynamics. We can compute /i for both fixed points B to 
find that 

18^ 
25' 



(32) 



6r' = 5r" , 



(30) 



hence, (a/2/3, 1, 0) is stable for /? > and 
{-y/2/3, -1, 0) is stable for /3 < 0. 

The fixed point C which corresponds to vanishing 
derivative of the potential. Since we have assumed that 
/3 — X is a constant, in this case the potential is flat 
(A = 0) in some finite region of X, instead of a point we 
actually have a curve C. Then its points live in an effec- 
tively two dimensional manifold, and one of the Lyaponov 
exponents is expected to vanish. Therefore we can now 
infer the nature of the critical point from the linear anal- 
ysis. Since the two nonzero eigenvalues are negative, we 
have what is called a local sink. Note that in fact A = 
corresponds to the massless field, and the reduction of 
the dimension of the phase space refiects the disappear- 
ance of a degree of freedom in the massless case due to 
the restored gauge invariance. 



7 



2. V ^ exp(-/3X2) 

This potential has a X"^ dependence hence is of the 
type we are looking for. If (3 is positive, it presents a 
maximum at X = and conversely, if /? is negative, it 
has a global minimum at the same value. 

Now the value of /i for fixed points B are 



(2) 



72 



(33) 



hence they are both stable if /3 > 0. 

The fixed point C is in this case {x,y,w) = (0,0,0) 
and its corresponding eigenvalues are mi^2 = ~(3/2)(l± 
\/l + 8/3/3) and 7713 = —37/2. This point is stable pro- 
vided /3 < 0, which clearly makes sense as it is the case 
that leads to a minimum in the potential. 



3. V = X^ + k 

Here A; is a positive constant introduced for the purpose 
of regularization purposes about x = 0. The quadratic 
potential is in a sense similar to the previous one with 
/3 < 0. Indeed we find that for fixed points B, 



(2) 



72 



1 



24 2/3 + fc ' 



(34) 



which mean that they are unstable. In this case, the 
fixed point C, also corresponds to {x,y,w) — (0,0,0) 
with eigenvalues 



TOl,2 



1T\ 1 



3fc 



1^3 = TT- 



thus this is a stable fixed point. 



4. V ^X* 



(35) 
(36) 



k 



Though the quartic potential seems very similar to the 
quadratic potential there are in fact some differences. 
Again k is a positive constant. For fixed point B wc 
compute 



(2) 



96 



1 



25 4/9 + ' 



(37) 



and again they are unstable. Fixed point C, which is 
again at {x,y,w) — (0,0,0) and also has eigenvalues 
(—3,0,-37/2) like fixed point B where the direction of 
vanishing eigenvalue is still given hy 6r — y^2/3 6x + Sy. 
Going to second order in perturbations along this direc- 
tion, like we did for point B, we find that /x^' =0. We 
thus have to go to the third order. 



(3) 



72 

125 /c ' 



(38) 



and there we find that the sign of the eigenvalue is neg- 
ative, thus, the point is a stable attractor. 



5. 1/=(X2-C2) +fc 

This potential has two minima ai X — ±C and a max- 
imum a.t X — (we are taking C > 0). For fixed point 
B we can readily calculate 



(2) 
/^B± 



144 (2/3 -C2) 
(2/3-C2)2 + fc 



hence, both fixed points are stable provided C > 
We have now, however, three type C fixed points: 



Ci: \^±C,± 

C2: (0,0,0). 
For Ci wc find eigenvalues 




mi^2 = I 1 ± 



24C2 



C,0 I , (40) 
(41) 

(42) 
(43) 

Thus, this fixed point is stable for k > 240'^ {C'^ - 2/Zf. 
We must also point out that this fixed point only exists 
for C < 1/2/3 as we must require \y\ < 1. For Ci we can 
compute the following eigenvalues: 



2 



ms = — — 



37 
2 



16C2 



3(C4 + k) 



TO3 



37 
2 



(44) 
(45) 



consequently we can conclude that this fixed point is un- 
stable which is not surprising given that it corresponds 
to the local maximum of the potential. 

The properties of the fixed points for all these forms of 
the potential are summarized in Table [TTl 

We are going to focus on two forms of the potential that 
will suffice to describe the general properties of this sys- 
tem that we have formally described above. In Fig. [5] we 
show the phase space portrait for a potential of the form 
V{X) — X^. It shows that the trajectories approach the 
minimum of the potential at X = and oscillate around 
this value. In Figs. [6] and [7] we show the phase space 
trajectories for the potential V{X) = {X^ — C^)^, where 
C is a constant. In the example of Fig. [51 C = 0.5 and 
we see that the late time behavior consists of the tra- 
jectories oscillating around the fixed point X — C and 
y = 0. Finally for the example of Fig. [71 we see that 
the trajectories approach the critical points x — ±-\/2/3, 
y = ±l. 
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V{X) 


A 


B 


C 


exp(-/3X) 


U 


B+ stable for /3 > 
B- stable for /3 < 


S 


exp(-/?X2) 


U 


stable for /3 > 


stable for < 


X' 


u 


U 


S 


X* 


u 


U 


s 


{X^-C^f 


u 


stable for C > ^2/3 


Ci stable, C2 unstable 



TABLE II: Stability of the fixed point in four classes of models. U - unstable; S - stable. 



V =X^ 



V ={X^ -C'^f, C = 15 




tarh(a;) 



FIG. 5: The phase space trajectories for quadratic potential. 
The minimum of the potential at X = has been marked 
with a dot. 




taili(a;) 



FIG. 7: The phase space trajectories V — (X^ — C^)^, when 
C = I > The stable fixed points at the minima of the 

effective potential have been marked with a dot. 



C=Q5 




tarh(a;) 



FIG. 6: The phase space trajectories V — {X^ — C^)'^, when 
C = ^ < \/2/3. The stable fixed points at the minima of 
the potential (both the true and the effective potential) have 
been marked with a dots. 



IV. COSMOLOGICAL PERTURBATIONS 

The general perturbations about the FLRW back- 
ground can be parameterized by writing the line element 
as 

+ a^{t){l-24))dx'dx^ + a^{t)hijx^x\ (46) 

where the two scalar perturbations -0 and <f> are the usual 
Bardeen potentials in the longitudinal gauge, hi is a 
transverse vector and /i.y is transverse and traceless as 
it describes the tensorial perturbations. 

The field equations for the scalar perturbations are 
then 



V2 



a 



V2 
— < 



'SH{^ + HxP) = -AnGSp, (47) 

(48) 



= 47rG'(p + p)- 
a 



0+ir(3(/.+V^) + (2i7+3i/2)^_ J_(0_y,) ^ 47rG(5p, (49) 



V_2 

3a2 
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^lJ) = 12ttG{p + p)w. 



(50) 



The first equation is the energy constraint (G'^q compo- 
nent), the second is the momentum constraint {G^^ com- 
ponent) involving the velocity perturbation 9, the third 
is the trace of the spatial part (G\ component), and the 
last one gives the shear propagation equations part (G'^ 
component) for the shear w. 



component, corresponding to the field itself, while the 
potential depends on the comoving field X. This variable 
K is just the dimensional version of the y we used in the 
phase space analysis. We notice also that the constraint 
equation ([4]), which may be written as 



Ft [ X 



Vx{aa + 3Xa)-Vx, 



(62) 



A. The three-form 

To parameterize the fluctuations of the three-form, we 
employ a similar decomposition as for the metric. The 
four degrees of freedom in a three-form then turn out to 
be two scalar and two vector degrees of freedom. The 
components of the three-form are fully specified by 



Aijk = a^{t)t^jk{X{t) + ao). 



(51) 
(52) 
(53) 



Here ak is a transverse vector and thus has two indepen- 
dent degrees of freedom. One easily sees that as usually, 
the vector and scalar perturbations decouple at linear 
order. The square of the field is then 



^2 ^ 6 [X^ + 2X{aQ + iXci 



(54) 



Under general gauge transformation x'^ ^ -f , spec- 
ified by the vector = (C°,C,i + f ), where f = 0, the 
field fluctuations transform as 



a — > a — aX£^, 
at ^ ai~ aX£^i . 



(55) 
(56) 
(57) 



Here we use the variable K = X + 3HX. In this trans- 
formation the metric potentials transform as 



(58) 
(59) 



The equations of motion for the scalar perturbations 
then are 



do + 3Hao 



( 

X 



-a + k{3(l)~ij)^0, (60) 



is not independent but follows consistently from Eqs. ([50]) 
and ([6T|) . The density, pressure and velocity perturba- 
tions come out as follows: 



Spx = K 



do + 3Hao 5-a + -ftr(30 — tp) 



+ Vx{ao + 'SX(j)), 



(63) 



Spx = -6px + iVxxX + Vx)iao + 3X(b), (64) 



{px +px)dx = jVxa. 



wx = 0. 



(65) 



(66) 



The first equation is the energy piece of the energy mo- 
mentum tensor Eq. ^ (T'q component), the second is 
the trace of the spatial part {T\ component), the third is 
obtained from the momentum containing part (T° com- 
ponent), and the last one gives the symmetric traceless 
part of the spatial T ^ components. Thus we get that the 
anisotropic stress due to the three-form vanishes. 

In the absence of other vector sources, the rotational 
perturbations evolve like 



h + Hh = 0, 
V2. 



V.x{Xh-a{). 



(67) 
(68) 



Thus the vector perturbations decay and can be ignored. 
Since the three-form does not generate tensor perturba- 
tions, their evolution equation is 



h + 3Hh ^/i = to(*\ 



(69) 



where ro^*) is tensorial anisotropic stress to which the 
three-form does not contribute. 



do + 'SHao + {3H + Vxx)ao- ^{a-^Ha) (61) 

+ k{?,4> - ^) + ?,{VxxX - Vx)(l3 + Wxi^ = 0. 

which represent the two independent components (Oij 
and ijk) of Eq. ([3]). The background equations of motion 
was used to simplify the second one. As discussed in 
the previous section, the kinetic terms are given by this 



B. Three- form domination 

Assume that the three-form dominates. Since it does 
not generate anisotropic stress, Eq. ([50]) tells us that 
-tp = (j). We can now derive an evolution equation for 
the Bardeen potential in a closed form. Eq. ([5T|l can 
be used to eliminate do from the system. Equations (j47p 
and (|63p may then be used to eliminate ao- Note that 
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only a appears always without derivatives in these equa- 
tions. We may eliminate a using Eqs. and ([55)1 . Fi- 
nally, plugging the solutions into Eq. (j49| with the right 
hand side given by (|64p we get 



propagation is a possibility. For positive however, the 
value of X is constrained to be < 1/2(3 to avoid an 
imaginary speed of sound and consequently an unstable 
scenario. 



H X 



(70) 



We can verify this result by differentiating Eqs. (|48|) and 
(j60p solving for a and ao, and verifying that the latter 
agrees with Eq. The RHS of Eq. ([TDD is simply 

6px /2 in the comoving gauge, and the density pertur- 
bation in the comoving gauge is given directly by the 
Poisson equation. Thus we read off the rest frame sound 
speed of the three-form: 



_ H X ^ 
' H X 



VxxX 
V 



X 



(71) 



where we have used the background relations in the three- 
form dominated universe. For a power law potential 
V{X) — X", it results that the speed of sound is a con- 
stant given by c| = n — 1. The expression (j7ip can also 
be found by noting that in the rest frame a = the 
expressions (|63 |) .(|64 p assume the form 



Spx\a=o 



1 



■Vx5A\ 



YIX 

^PxU=o = ^VxxSA"^. 



(72) 
(73) 



Thus the relation of the dynamical sound speed, Cg = 
6px I^Px J to the derivatives of the potential is general 
and not restricted to three-form dominated background. 

To analyze the behavior of the sound speed in more 
detail we have considered specific forms of the potential. 



Power-law potentials, V = X" 



3. Ginzburg-Landau potentials, V — {X"^ — C^)^ 

Now the expression for the sound speed reads 

2^2 



4 



1 



X2 - C2 ■ 



(74) 



In this case we find that the speed of sound is positive 
provided that either X"^ > or X^ < C^/3. Recall from 
the Section [mi that in the case C < ■y/2/3, the minimum 
of the potential a,t X — ±C is the late time attractor 
and that at this point the field changes its nature from 
phantom to non-phantom, or vice versa. At the level of 
the background kinematics, this is perfectly legitimate 
and can be easily verified by explicit solutions. However, 
scrutiny of the perturbation dynamics reveals that at the 
phantom divide the sound speed squared diverges and 
jumps from negative to positive infinity or the other way 
around. Clearly the phantom divide crossing is inhibited 
in reality for this particular potential. 

The form of the potential is such that it suggest that 
the origin corresponds to an unstable cosmological con- 
stant. This motivates us to consider a case where the 
three-form drives an inflationary period at the local max- 
imum of the potential at X w which can be seen as a 
natural initial value. As the field eventually drops to 
the true minimum, several possibilities may take place. 
If C > ^/2/3 then X is constrained to be \X\ < yj^ji 
hence for > we must require C > \/2, and the 
field will asymptotically approach X ~ ±-y/2/3. In the 
C < ^2/3 case, reaching the minimum will lead to insta- 
bilities which might be beneficial by generating efficient 
reheating. This possibility might be worth being ana- 
lyzed in the future. 



If the potential is = X", the sound speed squared is 
a constant c| = n — 1. The quadratic potential resembles 
a canonical scalar field as its speed of sound is always the 
speed of light. With higher order self- interaction terms 
the fluctuations in the three-form field propagate faster 
than light, which might be seen as a problem. Less con- 
troversial is the fact that potentials with n < 1 are un- 
stable and thus are not viable models of inflation. 

H. Exponential potentials, V — exp (— /9X^) 

Since the potential V — exp {—/3X'^) can be approx- 
imated with a quadratic correction to a constant term, 
the behavior of the sound speed is similar to the power- 
law case, Cg = 1 — 2f3X^. For negative f3, superluminal 



C. Scalar and tensor power spectra from inflation 

As we have seen, it turns out we can describe the scalar 
fluctuations of the field with only one degree of freedom 
by exploiting the constraints of the system, in particular 
Eq. dMl). This is due to the symmetries of the FLRW 
metric. One cannot see directly from the action that one 
of the four degrees of freedom present can be eliminated, 
as one should see when it is a general property of the 
theory. In Bianchi backgrounds, we have many examples 
of form solutions that are constrained. Of course, cosmo- 
logical perturbations in principle allow all the degrees of 
freedom present to propagate, but it happens now that 
the kinetic term has the gauge symmetry which reduces 
the number of physical degrees of freedom in the absence 
of the potential. Even when the potential is turned on. 
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the symmetry is partially efficient. This is because the 
potential depends only on and the spatial components 
of A are forced to vanish in the FLRW background so 
their fluctuations a cannot contribute at the linear order 
to the quadratic invariant A^ (see Ea. (|54|) ). Therefore 
we believe that the absence of the scalar second mode 
is due to linearization about the isotropic and homoge- 
neous solution. Higher order perturbations would thus 
be interesting to consider, but their study is beyond the 
present scope. 

To consider quantum fluctuations during inflation, we 
must find the canonical variable to describe the degree 
of freedom we have. It is conventional to refer to the 
curvature perturbation ( given by 



We must now solve the following equation of motion for 
the wave modes, 



wfe = 0. 



(81) 



We make the assumptions that the evolution of the uni- 
verse is power law like with scale factor a = {^t)p with 
p=— 1/(1 — e) which corresponds to 



(82) 



being approximately constant and that the evolution of 
the sound speed can be approximated with the power-law 
form, 



H- 



Hcj) 



that evolves as 



c 



H 



H , 1 



+ ( 



(75) 



(76) 



as can be verified using (|70p . Thus the curvature pertur- 
bation is conserved at large scales. By comparing with 
well-known cases in the literature 8l|, we can deduce 
that the canonical degree of freedom is now related to 
the curvature perturbation as 



-a^J-2H 
KCsH 



(77) 



At this point it is convenient to switch to conformal time, 
and in the remainder of this section a prime will denote 
derivative wrt to conformal time r = J adt. It is straight- 
forward to show that the canonical variable now obeys 
the equation of motion 



(78) 



The action for this variable v could be computed by ex- 
panding the action to second order, but this is not nec- 
essary. In the present case the equation of motion fixes 
the action, though only up to a constant. However, we 
know the normalization from analogy to previous litera- 
ture [sT]. Thus we may write 



I J ~ cl^''v,v, + '-^v^^ ^d^x, (79) 

where jij is the metric of spatial sections, which we here 
assume to be flat for simplicity. 

We proceed quantizing v by promoting the perturba- 
tion to an operator and expanding in plane waves. 



v{t, k) — 



d^k 
(27r)3/2 



t^fc(T)ake'^-^-f t^fe(T)a^ 



(80) 



(83) 



We will see shortly that this is in fact the case for the 
three-form during slow-roll (super)inflation of this work. 
In the case when we have z" /z cx t~^, the general solu- 
tion for u>k can be written as a sum of the Hankel func- 
tions, and the appropriately normalized solution with 
positive frequency in the asymptotic past is 



1 



a 



(84) 



where H'iF' is the Hankel function of the first kind of 
order z/, with 



^1 -4(-2"/0)t2 



1 



2(1 + a) 



(85) 
(86) 



Therefore the solution depend on the parameter cr, and 
presumably, trough z" / z^ on e and other slow roll param- 
eters. 

In the long wavelength limit, a; <C 1, the behavior 
of the Hankel function is h';^\x) {i/'K)T{v){2lxY . 
From this this asymptotic behavior we can calculate the 
power spectrum of the curvature perturbation on the 
large scales 



27r2 
k 



Ai 



k 



(1 



14 

2v-l 



"s — 1 



l-2i/ 2 
CS K 



a ) 327r2 
Y(vf 



r(3/2)2 



(87) 



(88) 



We can then read the scalar spectral index which is 

ns-l = 3-2i/. (89) 

Because the three-form does not introduce tensor 
sources, and the evolution of gravity waves is given by 
the usual equation (j69p . we obtain in the standard way 
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the spectrum of tensorial perturbations for power law 
evolution of the universe, 



Vt 



k 



3-2m ^2 



where 



27r2 

22^-1 (1 _ g)2A<- 



1 np^? „2 

r(3/2)2 ' 



7Tv/l-4(-a"/a)r2 



The tensor spectral index is 



(90) 



(91) 



(92) 



We will now compute these quantities in terms of the 
slow-roll parameter which can be calculated for a given 
scalar potential. 



1. Slow-roll (super)inflation 

In order to proceed with the computation of in 
Eq. we note that, 



2 a\ 



(93) 



where the slow-roll parameter e in our system is written 
as 

In what follows we mean inflation when e is positive and 
super-inflation when it takes negatives values. In order 
to obtain inflation we do not need slow-roll in the sense of 
negligible velocity of X (in fact a flat potential leads to de 
Sitter inflation even if the field is moving, see Eq. (I14p ). 
Here, however, we are now interested in investigating pre- 
cisely the slow-roll case. Then the velocity of X can be 
neglected, and thus we find that e can be well approxi- 
mated by 



3 _Kx 
2 V 



X 



(95) 



which allows us to immediately determine whether the 
universe is inflating for a given choice of the scalar po- 
tential at a given value oi X. Using the original equation 
of motion for X Eq. ^ and neglecting the X contribu- 
tion, it can be found that 



(96) 



This can also be written, using (|95p in the useful form 



V 



X 



4 



N 



9 (kX)2 



(97) 



Differentiating e in Eq. ([94]) with respect to conformal 
time and using the equation of motion for X it is obtained 

e'^2e^riaH, (98) 

where we have used Eq. ([57)1 and defined 



1 - - 



9 {nXy \ 



( Vxx 



V 



V 



y2 



VxX 



(99) 



We therefore see that the second term in 77 is suppressed 
by e and typically 77 is of order unity. However, Eq. ([98]) 
tells us that dlne/dN = 2eri and suggests that e can 
be considered constant when small or equivalently, that 
the evolution of the universe is power law like with scale 
factor a — {—t)p and p = —1/(1 — e). 

Similarly, by differentiating the speed of sound, it is 
obtained that 



c's — e 9 cs aH . 
where we have defined 



(100) 



9; 



1 



9 {nXf 



V,xxV 

,X 



1 - 



V 



XXX 



^ .XX 



Vx \ 
VxxX j 



(101) 

For a power law potential we find that 9 vanishes and 
therefore cs is a constant in agreement with what we 
found earlier. Therefore, when e is small, both e and 9 
are nearly constant which enables us to indeed write (j83[) 
with 



Now Eq. ([89| gives that 



ns-\ 



-2e ( 1 + -77 



(102) 



(103) 



We immediately see that near and approaching the fixed 
point B where e is negative, the three-form predicts a 
slightly blue spectral index for curvature perturbations 
which is disfavored by current observations. This is the 
case illustrated in Fig. 21 When e is positive, however, 
a red tilted spectrum is obtained. Such an example is 
shown in Fig. [3l 

From Eq. ([5^ the tensor spectral index is equivalent 
to standard scalar field inflation 



TiT = — 2e . 



(104) 



Like the curvature spectral index, this is predicted to be 
slightly blue when the evolution is near and approaching 
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the fixed point B, and red tilted otherwise. The tensor 
to scalar ratio is, however, modified. For small e we have 



4 



16cs|e| 



(105) 



Thus it is in principle possible to distinguish the three- 
form inflation from scalar field already from the spectra 
of linear perturbations. 



D. Matter perturbations 

The impact of the presence of three-from to the evo- 
lution of matter inhomogeneities is considered here. The 
observationally relevant case is dust-like matter (cold 
dark matter and baryons) in the late-time universe, 
where anisotropic stresses can be neglected. The con- 
servation equations for matter then are, in the Fourier 
space, 



3*, 



a 



= —H9„ 



fc2 



Combining these gives 



2H5„^ + —(j)^3(l) + 6H(j) . 



(106) 
(107) 

(108) 



The effect of the three-form comes thus through the back- 
ground expansion (H and a) and through the coupling 
of the matter to the gravitational potentials. We then 
specialize to the subhorizon scales. Then the RHS of 
the previous equation can be neglected. We need yet to 
evaluate the gradient of the gravitational potential. At 
this small-scale limit the perturbed energy constraint is, 
using (gT]) and 



1 



Pin ^771 

+ Vx l^a + ao + 5xA . (109) 



By using the momentum constraint (|48p with (|S5)) and 
d07\ . we get 



Vxa = 2{(b + H4>) + -pUSm - 3?i) . 



(110) 



This shows that the contribution from a to the gravi- 
tational potential in (|109p is suppressed by the a^/k'^, 
which allows us to neglect it at the small scale limit, at 
least excluding the special case of diverging K/X. As- 
suming {k/a)^ ^ V_x/X, we can furthermore show that 
the equations of motion (|60p and (|6ip imply, instead of 
([62)) . the algebraic relation 



between the scalar perturbation and the gravitational po- 
tential. Using this in (|109p and plugging back to (jlOSp . 
we get 



2H5„, = ATiG{k,t)pmS,, 



(112) 



where the scale and time dependent effective gravita- 
tional constant is given by 



8TrG{k,t) = + 



1 VxXa^ 

2 c|/c2 • 



(113) 



It is to be expected that the apparent modification is 
proportional to the slope of the potential, since in this 
case this slope vanishes the field reduces to a smooth cos- 
mological constant. We also immediately note that if the 
sound speed squared is large, the small scale limit corre- 
sponds to the presence of a smooth component. Possible 
new effects (other than due to modified background ex- 
pansion) require sound speed less than unity, and the rel- 
evant range of scales is the intermediate regime between 
the cosmological horizon and the sound horizon: 



k' 



(114) 



Vxxaa = -{iV.xxX + V.x) 



(111) 



The possible scale-dependent signature extends also to 
larger scales, but from superhorizon scales one typically 
expects possible observable effects only to the largest cos- 
mic variance limited CMB multipolcs. For small sound 
speeds the modifications can extend to nonlinear scales, 
where we cannot trust the expression ()113p anymore. In 
the regime we can trust it, it is in principle possible to 
constrain some classes of three-form cosmologies by us- 
ing the probes of large scale structure, weak lensing, the 
integrated Sachs- Wolfe effect and their correlations. 



V. FORMALITIES 

In the following two subsections we will consider some 
formal aspects of our simple theory. The reader inter- 
ested solely in cosmological phenomenology may proceed 
to the conclusions. In particular, we will scrutize the 
field content of the action ^ in terms of a so called 
Stiickelberg trick [s^ and by performing mappings and 
duality transformations to other theories. One of the 
aims is to proof explicitly the claim made in the introduc- 
tion that the models can be seen as a novel description of 
scalar field models in some specific situations but how- 
ever for most cases it is not (even formally) equivalent 
to a scalar theory. This point, as well as some clarifi- 
cations of our terminology will be useful to make here, 
since many different conventions (and some confusions) 
exist in the literature. 



A. Gauge invariance and stability 

The antisymmetrized gradient of the three-form gauge 
potential term in Eq. ([1]) is gauge-invariant under trans- 
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formation A ^ A + [V A] , where A is an arbitrary two- 
form. However, the potential in the action ([1]) breaks this 
symmetry. This can be seen to result in extra degrees of 
freedom in the model, analogous to the appearance of 
longitudinal polarization of the massive photon in the 
Proca theory. To make this explicit, one may introduce 
a Stiickelberg form S in such a way that for a redefined 
A = A + 4[\/T,], the Lagrangian 



grangian l|119p 



C 



V{iA + F{T.)y 



(115) 



is then manifestly invariant under the gauge transforma- 
tions 



i^i+[VA], 



A/4. 



(116) 



In the S = gauge we recover our original Lagrangian. If 
there is a gauge where the gradient of A is orthogonal to 
A, and if we assume we can expand around a background 
solution given by constant A, we can write 

£ = £-V'{A^)F^{J:), (117) 

making transparent the appearance of the extra two- form 
degree of freedom. This also seems to imply that the ex- 
tra degree of freedom becomes a ghost when V'{A'^) < 0. 
We remind that this is also exactly the condition for the 
three-form to violate the null energy condition in FLRW 
background, i.e. to become phantom-like with equation 
of state less than —1. However, investigating the stability 
of the field in more detail by considering inhomogeneous 
and anisotropic fluctuations of the field and taking into 
account their backreaction due to coupling to the metric, 
we find that the conditions for the stability of the canoni- 
cal degrees of freedom can be more subtle than the naive 
implication of (|117p . In fact, the behavior of the field 
depends on the second derivative of the potential since 
the propagation speed of the physical fluctuations turned 
out to be given by (equivalently to the expression (|7ip in 
terms of the comoving field) 

The field, even if phantom, can be stable at least classi- 
cally. However, divergences tend to occur at the "phan- 
tom divide" when V'{A'^) crosses zero. Similar phe- 
nomena, linking classical singularities and quantum no- 
ghost conditions have been observed in other models 
[27l . [ssl . [s^ . The conditions for the possibility of a viable 
phantom crossing is now the following: twice differen- 
tiable V{x) exists for positive x in such a way that V'{x) 
changes sign at x = xq, and xoOj^{V' (xq)))' is finite. 

This can easily change if nonminimal couplings arc in- 
troduced, but in the present paper we consider only min- 
imal couplings. 



B. Duality and equivalence of theories 

In this subsection we will discuss some equivalences 
between n-form models. First we introduce a parent La- 
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Cp = -F'--AV-F-ViA^), (119) 
which can be rewritten in terms of a Faraday form 



Cf=f{F\x))~V{x')., 



(120) 



where / and V are arbitrary functions and x is a n-form, 
which describes our starting point JT]) for a: = A. We 
will also show that also a four-form can emerge from it, 
i.e., we can rewrite the parent Lagrangian in terms of 
gauge-fixing terms only 



^9 = 9 ((V 



U{x^), 



(121) 



where g and U are functions. Next we show how the 
Hodge dual of the parent Lagrangian (|126p can give rise 
to a vector or a scalar field description. Finally, the chain 
of equivalences is summarized in the diagram ([8]). Du- 
alities in the case of nonminimal gravity couplings have 
been discussed in [s^]- 

Starting with the parent Lagrangian pi9p . where F 
is an independent four-form and solving its equation of 
motion, one gets F = —4[VA]. By plugging this back 
into pi9p one obtains the original action ([T]). 

One may also integrate out the three-form and obtain 
a theory for the four-form F as follows. Varying with 
respect to A gives us the equation of motion 



V ■F + 2V'{A^)A = 0, 



implying 



y = (V • F)2 ^ 144 {V'{A^)y A^ 



(122) 



(123) 



If we now plug A from and the solution A^{Y) 
into (|119p . we get a dynamical four- form theory: 



48 



-F' 



Y 



72 V (A2(r)) 



V[A^{J)). (124) 



The Lagrangian is now written in a gauge-fixing form 
(|12ip with X = F. The Faraday form constructed from 
a four-form is of course trivial, F{F) — 0, and static 
four-forms contribute only a constant. As mentioned in 
the introduction, this has been employed in attempts to 
solve the cosmological constant problem [t^ [t^ . Four- 
form formulation of f{R) gravity [s^l and dark energy 
from promoting the Levi-Civita symbol into dynamical 
form have been considered recently [i^j. 

The Hodge dual of the three-form is a vector {*A). 
Writing the parent Lagrangian (|119p in terms of the dual 
forms 

= e(*F) EE e$ , A = e{'fcA) = eB, (125) 

where $ is a scalar field and i? is a vector, we obtain 

£ = -i$2 _ 5 V$ _ y . (126) 
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Now the equation of motion for $ is simply that <I> = 
V • i3, and replacing this back gives us the self-coupled 
vector theory 



£ = i(V-B)^-y(-BV6) 



(127) 



in a gauge-fixing form as in (|12ip with x = B. Re- 
cently the cosmological significance of the Maxwell the- 
ory supplemented with the gauge- fixing term like in (jl27p 
has been considered, and very interestingly it has been 
found that the gauge-fixing term results in an effective 
cosmological constant in the curved background (or al- 
most constant, since it fluctuates) while the theory in the 
Minkowski limit reduces to standard electromagnetism 
[ssl . l86l |. Such a Maxwell theory is dual to three- form 
having both the F'^{A) and (V • A)^, but presently we 
confine to include only the canonical term yielding only 
the gauge- fixing term for the corresponding vector. 

Finally, we may integrate out the vector from (jl26p to 
obtain a scalar field theory. The Euler-Lagrange equation 
B is 



implying 



V$ = -V^' {-B^/6)B, 



9{V<i>f^[V' {-B^/6)]'b' 



(128) 



(129) 



Similarly as with the four-form, we assume these equa- 
tions are invertible, and then write the Lagrangian (I126P 
in terms of them as 



2 3V'i~B^Ys)/6) 



-V{-B\Ys)/6), (130) 



where now Ys = (V<i>)^. The dual Lagrangian is now 
written in a Faraday type (|120p with a; = This com- 
pletes our task of deriving the equivalent reformulations 
we mentioned in the introduction. 

It is useful to note that only two possibilities of canon- 
ical forms exist in four dimensions. Any such form is 
either a vector or a dual vector. The former is dual to 
a three-form, and can be rewritten as two-form which is 
self-dual. The latter can be seen as scalar field, thus dual 
to a four-form and consequently rewritable as a three- 
form. Schematically, the chains of equivalences can be 
written as 



/I 



53 = 

/o 



> gA 



g2- 
/3 



/I, 
>5l- 



(131) 
(132) 



Here /n denotes an n-form described by a Faraday type 
Lagrangian (|120p , and gn an n-form described by gauge- 
fixing terms type Lagrangian (|12ip . Duality in the Hodge 
sense is indicated with <-!■, and = means equivalence be- 
tween Lagrangians by change of variables. This is also 
shown in diagram [S] 

The case we consider in the present paper belongs to 
the chain p32p or the third group in diagram [51 Thus it 
seems we are, at the same time, discussing four types of 
theories: 



g 





B = B(V«C) . 
- > 



C = C(VB) 



^^C=_C(V.A)_^ 
> 



A = A(V«F) 4 
- >■ 



4 
f 



3 

g 



A = A(VC) 



F = F(VA) 



* = <I>(V«B) 

> 

B = B(V<I>) Q 
f 



1 

g 



FIG. 8: A summary of the chains of equivalences between 
forms. In the case of canonic form, / refers to a model with 
the Faraday kinetic term with an action of the form (|120p and 
a; to a model with the dual (gauge-fixing) kinetic term with 
an action of the form H12ip . The Hodge duality operation 
is vertical movement denoted with a star, and a horizontal 
movement is a change variables for which schematic formu- 
las are given. The first group in the figure, consisting of the 
form-form and dual, is trivial. The second group, which de- 
picts the chain (|13ip . would correspond to a starting point of 
having the gauge-fixing kinetic term for the three-form. The 
third group depicts the chain (|132p . corresponding to the case 
whose generalization we consider in this paper. 



• a model of vector with self-interactions and a 
gauge- fixing type term, Eq. (|127p . 

• a K-essence type scalar field model Eq. (|130p . 

• a four-form with dynamics due to a nonstandard 
kinetic term, Eq. (|124p . 

• a canonical three-form, Eq.((T]). 

These equivalences are valid due to the possibility of 
rewriting the four degrees of freedom in the three-form as 
a vector. Furthermore, if the vector is exact, i.e. express- 
ible as a gradient of a scalar, the model can be reduced 
to a scalar field, equivalently a four-form. 

However, we should immediately mention that these 
equivalences break down in many cases. In the above 
derivations, this breakdown occurs when there are no 
real solutions to equations (|123p or (|129p . The exam- 
ple potential of most interest in the present study of 
the action ([T]) is of the form of a displaced power-law 
ViA^) = {A^/G - C^)^, since that is the case that in- 
cludes most of the features of the possible dynamics. In 
particular, having a non-nionotonic first derivative allows 
the field to dynamically change its nature. A scalar field 
formulation, however, seems not to be available. An- 
other simple class of potentials we considered is the ex- 
ponential potential. As shown in appendix [XI this class 
does not seem to admit a scalar field formulation either. 
Even in some cases where the equivalence formally holds, 
the formulation in the three-form language is much more 
transparent. 
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C. Generalizations 

One might consider extensions of the simple and min- 
imal action ([T]). Since, due to the self- interactions of the 
gauge potential A, there is no gauge invariance, the gauge 
invariance cannot be used as a criterion for the action. 
However, gauge-invariant nonlinear terms can be consid- 
ered and motivated by quantum corrections such as in 
nonlinear electrodynamics [13, HI] ■ Here a phenomeno- 
logical motivation for such terms was briefly mentioned 
at the end of section[ITl When going to curved spacetime, 
it is known that terms involving contractions of F and the 
Riemann and Ricci tensors arise as leading order quan- 
tum corrections [s^. Models obtained by the inclusion 
of the gauge-invariance preserving forms like f{R, G)F^, 
where f{R, G) is some function of the curvature scalar or 
Gauss-Bonnet term, have been applied in cosmology in 
cases of Maxwell or Yang-Mills fields [13, [sH, HI • 

As mentioned, the potentials break gauge-invariance 
and one can allow more gauge-invariance breaking terms. 
These could still be quadratic. In particular, a gen- 
eral quadratic second order Lagrangian for a three-form 
would read, 

+ a^A- Ric- A + a^A- -{n- -A), (133) 

where TZ is the Riemann and Ric the Ricci tensor. There 
are thus six coefficients to specify such a theory. If the 
principle of minimal coupling is kept as a guide as in the 
present study, only three are left. A particular combi- 
nation of the coefficients, keeping only a2 = but fixed 
values for the rest ai results in a scalar-field like equation 
of motion for the comoving field in a FLRW background. 
For this reason such a model has been considered as gen- 
eralization of the vector inflation scheme [13, Hli- In 
future work, it would be interesting to study the con- 
sistency constrains of nonzero a^, in particular the im- 
plications of the gauge-fixing term given by a nonzero 
a2- Many of the parameter possible combinations in the 
general action (|133p are probably excluded due to ap- 
pearance of ghosts and instabilities. In addition. Solar 
system experiments can be used to constrain the non- 
minimal couplings to gravity, as, Q!4 and a^. To our 
knowledge these issues have not been addressed in the 
case of three-form, whereas vector models have been ex- 
tensively studied (see Ref. ^] and references mentioned 
in the introduction section Finally, since the Maxwell 
field allows a dual description as a three-form, one may 
contemplate whether it is possible that nonminimal cou- 
plings can lead to variations of fundamental parameters 
such as the fine structure constant. 

VI. CONCLUSIONS 

We considered the evolution of the universe in the pres- 
ence of three-forms. We assumed a canonical and min- 



imally coupled action taking into account the possibil- 
ity of self-interactions of the form field. Then a form 
with three differing spatial indices is compatible with an 
isotropic and homogeneous cosmological background. It 
turns out that such form, despite its canonical form, quite 
generically violates the usual energy conditions. The 
strong energy condition is violated when 

VxX < ^ Q(i: + 3HXf + ^ wx < 

. 

This happens quite easily. Slow roll is not required for 
accelerating behavior, only that VxX is not large com- 
pared to the energy density. The null energy condition is 
broken and the field behaves as phantom when 

V,xX < => wx <-l. (135) 

Generally stability problems appear in crossing the phan- 
tom divide, but they might be overcome by nonminimal 
couplings. In the absence of a potential, the effect of the 
field reduces simply to generating a cosmological con- 
stant. Thus the three- form seems a very suitable culprit 
for the accelerating expansion which is believed to take 
place both at an early stage and at a late stage of the 
history of the universe. One notes that phantom divide 
crossing is possible for simple forms of the potential, such 
as the Landau-Ginzburg form that we took as our main 
example. 

We performed a phase space analysis of the model and 
found three distinct fixed points: matter domination, 
a potential extremum and a peculiar fixed point corre- 
sponding to a kinetic domination of the three-form (and 
a potential domination of the dual scalar field Eq. (|130p 
in the special cases that such a scalar field exists). The 
latter two fixed points describe de Sitter spaces, and their 
nature and stability depends on the potential. The math- 
ematical properties of the system proved nontrivial, and 
forced us to go up to third order in the perturbations for 
some forms of the potential. This tempts one to investi- 
gate the model in the bifurcation theory framework. On 
the physical side, one of the de Sitter fixed points is al- 
ways an attractor. Thus the fixed points are interesting 
for inflationary and dark energy applications. For the 
latter, one wishes to find initial conditions independent 
dynamics before the acceleration in order to reduce the 
fine-tuning. Previously it was shown that scaling and 
tracking solutions exits and the requirements for these 
to last for a large number of e-folds can be quantified 
[T^I- It is then clear that one may realize a rich variety 
of background dynamics using a single three-form as the 
energy source. 

The three-form fluctuations were then investigated. 
They were parameterized in terms of two scalar and two 
vector modes (whose form was motivated by the dual 
vector). The latter are (as usual) phenomenologically 
less interesting. However, it is worthwhile to note that 
in general these degrees of freedom also exist. It turned 
out that one may describe the scalar perturbations with 
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one effective degree of freedom, since the constraints al- 
lowed to eliminate one of the two fields parameterizing 
the scalar fluctuations. The presence of the three- form 
can modify the evolution of scalar fluctuations in matter 
distribution. If the three-form sound speed is sufRciently 
less than unity, one expects possibly observable effects 
for the large scale structure, weak lensing, the integrated 
Sachs- Wolfe effect and their correlations. The impact of 
three-form fluctuations can be quantifled by introducing 
an effective gravitational constant defined by 



8TTG{k,t) 



2 Vxxk^ 



(136) 



During a three-form driven inflation the nontrivial dy- 
namics can also lead to sound speed dependent possibly 
detectable signatures. We identified the canonical de- 
grees of freedom and quantized them. A slow-roll param- 
eterization was reconsidered for this new case and a very 
convenient way of studying the quantum generation of 
perturbations near a de Sitter fixed point was developed. 
The spectral indices of scalar and tensor perturbations 
have a easily computable dependence on the form of the 
potential through the sound speed. 



VxxX 
Vx 



(137) 



The tensor to scalar ratio is modified directly by this 
quantity. 

The four formulas (|134p ~ (|137p summarize how the 
shape of the potential determines the nature of the field 
and thus the background dynamics and properties of fluc- 
tuations for a given model. To conclude, we have shown 
that these objects, also present in string theory, can give 
rise to viable cosmological scenarios with potentially ob- 
servable signatures distinct from standard single scalar 
field models. 



APPENDIX A: THE SCALAR FIELD 
FORMULATION 

In this Appendix we derive explicitly the scalar field 
formulation corresponding to a simple exponential po- 
tential and a displaced power-law potential. The main 
conclusion to be drawn is that the scalar field formula- 
tion of these models is typically non-trivial. The map- 
ping of these models is not a bijection, but the defining 
equations have multiple solutions and branches. Further- 
more, some of the solutions become complex, so mapping 
becomes ill-defined or ceases to exist. 

Consider the potential V{A^) = Vq{<oA^ - C^f, where 
we have introduced the irrelevant rescaling by the factor 
of 6 just simplify some formulas. Then Eq. (|129p reads 



This a third order equation for the square of the form: 
in general three solutions exist. For the sake of explic- 
itness, we pick up the solution which seems to lead to 
the simplest form of the Lagrangian. This Lagrangian is 
then 



3K)(V$)2 



^/((V*)^) 



1^2 
^0 



where 
/((V$)2) 



32CV^ + 27V^{V^f 



-m2$2, (A2) 



(A3) 



3\/3 JVo**(V$)2 (27(V$)2 - MCWg 



This is a noncanonical field indeed. 

Consider the exponential potential V{A-^) = V^e'^^^ 
Then Eq.((l29l) reads 



(A4) 



The solution can be written in terms of the Lambert's 
VF-function as 



A' 



1 

2A 



where we defined 



-(V$)^ 
2\Vi 



(A5) 



(A6) 



Lambert's W-function has multiple branches. It is real 
W{x) € M, when the argument x> 1/e. The scalar field 
Lagrangian would then read 



(A7) 



It seems that in general such a model is not well defined. 
In any case it is clear that the formulation as a canonic 
form with an exponential potential is considerably more 
tractable than this formulation. 



APPENDIX B: ABOUT THE DUAL VECTOR 
FIELD 

Here we state some facts concerning the dual vector 
field *A appearing in (|125p . We obtain it by the change 
of variables, and it can be written explicitly using the 
indices as 



1 



(V$)2 = 4Vo2(6A2 - C'^fA^ 



(Al) 



*A„ 



(Bl) 



18 



Since the equations come out in equivalent form regard- 
less whether we begin with A or with *A, we find no ben- 
efit in switching between alternative descriptions. How- 
ever, one exception appears in our parameterization of 
the three-form perturbation in Eqs. (|5ip and ((52)) . 

Aoij = a^(t)eyfc(Q;,fe -I- Qffc), (B2) 
A,jk = a'^{t)e,jk{X{t) + ao). (B3) 

To justify calling the components of scalar perturbations 
as if they were components of a vector, one may employ 
the gauge transformations Eqs. ([55)) - ([Ff|) to see that they 
indeed transform as if they form a vector. A more direct 
way to see this is to note that the dual vector is just 

*Ao - + ao + 6X(l)), (B4) 

*Ai = ^ {a,i + ai) . (B5) 



One could have absorbed 6X0 into the definition of ao. 



ACKNOWLEDGMENTS 



We acknowledge n-formal discussions with Jose 
Beltran Jimenez, David Mota, David Mulryne, Sami 
Nurmi and Cyril Pitrou. The authors are supported by 
Deutsche Forschungsgemeinschaft, project TRR33. This 
work was initiated at the workshop "New horizons for 
modern cosmology" at the Galileo Galilei Institute in 
Florence; we would like to thank the institute for hos- 
pitality. 



G. Nordstrom, Phys. Z. 15, 504 (1914), physics/0702221. [23 
C. Brans and R. H. Dicke, Phys. Rev. 124, 925 (1961). 

A. A. Starobinsky, Phys. Lett. B91, 99 (1980). [24 
C. Wetterich, Nucl. Phys. B302, 668 (1988). 
P. J. E. Peebles and B. Ratra, Astrophys. J. 325, L17 [25 
(1988). 

L. Amendola, Phys. Rev. D62, 043511 (2000), astro- [26 
ph/9908023. 

T. Koivisto and D. F. Mota, Phys. Lett. B644, 104 [27 
(2007), astro-ph/0606078. [28^ 

B. A. Bassett, S. Tsujikawa, and D. Wands, Rev. Mod. 
Phys. 78, 537 (2006), astro-ph/0507632. 

E. J. Copeland, M. Sami, and S. Tsujikawa, Int. J. Mod. [29 
Phys. D15, 1753 (2006), hep-th/0603057. 
L. H. Ford, Phys. Rev. D40, 967 (1989). [30 
T. S. Koivisto and D. F. Mota, JCAP 0808, 021 (2008), 
0805.4229. [31 

A. Golovnev, V. Mukhanov, and V. Vanchurin, JCAP 
0806, 009 (2008), 0802.2068. [32 

B. Himmetoglu, C. R. Contaldi, and M. Peloso, Phys. 

Rev. Lett. 102, 111301 (2009), 0809.2779. [33 
M. Karciauskas, K. Dimopoulos, and D. H. Lyth (2008), [34 
0812.0264. 

A. Golovnev and V. Vanchurin, Phys. Rev. D79, 103524 [35 
(2009), 0903.2977. 

M.-a. Watanabe, S. Kanno, and J. Soda (2009), [36 
0902.2833. 

S. Yokoyama and J. Soda, JCAP 0808, 005 (2008), [37 
0805.4265. 

K. Dimopoulos, M. Karciauskas, D. H. Lyth, and Y. Ro- [38 
driguez, JCAP 0905, 013 (2009), 0809.1055. 
K. Dimopoulos, M. Karciauskas, and J. M. Wagstaff [39 
(2009), 0907.1838. 

V. V. Kiselev, Class. Quant. Grav. 21, 3323 (2004), gr- [40 
qc/0402095. 

C. Armendariz-Picon, JCAP 0407, 007 (2004), astro- [41 
ph/0405267. 

T. Koivisto and D. F. Mota, Phys. Rev. D73, 083502 [42 
(2006), astro-ph/0512135. [43 



C. G. Boehmer and T. Harko, Eur. Phys. J. C50, 423 
(2007), gr-qc/0701029. 

T. Koivisto and D. F. Mota, Astrophys. J. 679, 1 (2008), 
0707.0279. 

H. Wei and R.-G. Cai, Phys. Rev. D73, 083002 (2006), 
astro-ph/0603052. 

T. Koivisto and D. F. Mota, JCAP 0806, 018 (2008), 
0801.3676. 

J. B. Jimenez and A. L. Maroto (2009), 0905.1245. 

D. F. Mota, J. R. Kristiansen, T. Koivisto, and N. E. 
Groeneboom, Mon. Not. Roy. Astron. Soc. 382, 793 

(2007) , 0708.0830. 

T. R. Dulaney, M. I. Gresham, and M. B. Wise, Phys. 

Rev. D77, 083510 (2008), 0801.2950. 

C. Germani and A. Kehagias, JCAP 0903, 028 (2009), 

0902.3667. 

T. Kobayashi and S. Yokoyama, JCAP 0905, 004 (2009), 
0903.2769. 

T. S. Koivisto, D. F. Mota, and C. Pitrou (2009), 
0903.4158. 

B. Saha, Phys. Rev. D74, 124030 (2006). 

C. G. Boehmer and D. F. Mota, Phys. Lett. B663, 168 

(2008) , 0710.2003. 

W. Zhao and Y. Zhang, Class. Quant. Grav. 23, 3405 
(2006), astro-ph/0510356. 

K. Bamba, S. Nojiri, and S. D. Odintsov, Phys. Rev. 
D77, 123532 (2008), 0803.3384. 

E. J. Copeland, A. Lahiri, and D. Wands, Phys. Rev. 
D51, 1569 (1995), hep-th/9410136. 

A. Lukas, B. A. Ovrut, and D. Waldram, Nucl. Phys. 
B495, 365 (1997), hep-th/9610238. 

M. Gasperini and G. Veneziano, Phys. Rev. D59, 043503 
(1999), hep-ph/9806327. 

N. Bilic, G. B. Tupper, and R. D. Viollier (2005), hep- 
th/0504082. 

N. Bilic, G. B. Tupper, and R. D. Viollier, J. Phys. A40, 
6877 (2007), gr-qc/0610104. 

G. De Risi, Phys. Rev. D77, 044030 (2008), 0711.3781. 
J. W. Moffat, Phys. Lett. B355, 447 (1995), gr- 



19 



qc/9411006. 

[44] T. Damour, S. Deser, and J. G. McCarthy, Phys. Rev. 

D47, 1541 (1993), gr-qc/9207003. 
[45] T. Prokopec and W. Valkenburg, Phys. Lett. B636, 1 

(2006), astro-ph/0503289. 
[46] T. Janssen and T. Prokopec, Class. Quant. Grav. 23, 

4967 (2006), gr-qc/0604094. 
[47] P. D. Gupta (2009), 0905.1621. 

[48] E. I. Guendelman and A. B. Kaganovich, Class. Quant. 

Grav. 25, 235015 (2008), 0804.1278. 
[49] E. I. Guendelman and A. B. Kaganovich, Phys. Rev. 

D75, 083505 (2007), gr-qc/0607111. 
[50] E. I. Guendelman, Mod. Phys. Lett. A14, 1043 (1999), 

gr-qc/9901017. 

[51] E. I. Guendelman and O. Katz, Class. Quant. Grav. 20, 

1715 (2003), gr-qc/0211095. 
[52] C. Armendariz-Picon and V. Duvvuri, Class. Quant. 

Grav. 21, 2011 (2004), hep-th/0305237. 
[53] U. Gunther, P. Moniz, and A. Zhuk, Phys. Rev. D68, 

044010 (2003), hep-th/0303023. 
[54] T. Saidov and A. Zhuk, Phys. Rev. D75, 084037 (2007), 

hep-th/0612217. 
[55] T. Jacobson and D. Mattingly, Phys. Rev. D64, 024028 

(2001), gr-qc/0007031. 
[56] V. A. Kostelecky and R. Lehnert, Phys. Rev. D63, 

065008 (2001), hep-th/0012060. 
[57] B. Li, D. Fonseca Mota, and J. D. Barrow, Phys. Rev. 

D77, 024032 (2008), 0709.4581. 
[58] J. A. Zuntz, P. G. Ferreira, and T. G. Zlosnik, Phys. Rev. 

Lett. 101, 261102 (2008), 0808.1824. 
[59] S. M. Carroll, T. R. Dulaney, M. I. Gresham, and 

H. Tam, Phys. Rev. D79, 065012 (2009), 0812.1050. 
[60] L. Campanelh (2009), 0907.3703. 

[61] A. de Oliveira-Costa, M. Tegmark, M. Zaldarriaga, and 
A. Hamihon, Phys. Rev. D69, 063516 (2004), astro- 
ph/0307282. 

[62] D. J. Schwarz, G. D. Starkman, D. Huterer, and 
C. J. Copi, Phys. Rev. Lett. 93, 221301 (2004), astro- 
ph/0403353. 

[63] H. K. Eriksen, A. J. Banday, K. M. Gorski, and P. B. 

Lilje, Astrophys. J. 622, 58 (2005), astro-ph/0407271. 
[64] P. K. Samal, R. Saha, P. Jain, and J. P. Ralston, Mon. 

Not. Roy. Astron. Soc. 385, 1718 (2008), 0708.2816. 
[65] J. Hoftuft et al., Astrophys. J. 699, 985 (2009), 

0903.1229. 

[66] M. Blomqvist, E. MortseU, and S. Nobih, JCAP 0806, 

027 (2008), 0806.0496. 
[67] A. R. Cooray, D. E. Holz, and R. Caldwell (2008), 

0812.0376. 



[68] T. Koivisto, Phys. Rev. D76, 043527 (2007), 0706.0974. 
[69] C. Quercellini, M. Quartin, and L. Amendola, Phys. Rev. 

Lett. 102, 151302 (2009), 0809.3675. 
[70] C. Quercellini, P. Cabella, L. Amendola, M. Quartin, and 

A. Balbi (2009), 0905.4853. 
[71] M. J. Duff and P. van Nieuwenhuizen, Phys. Lett. B94, 

179 (1980). 

[72] S. W. Hawking, Phys. Lett. B134, 403 (1984). 

[73] N. Turok and S. W. Hawking, Phys. Lett. B432, 271 

(1998), hep-th/9803156. 
[74] T. S. Koivisto and N. J. Nunes (2009), 0907.3883. 
[75] S. S. Gubser (2000), hep-th/0010010. 
[76] A. R. Prey and A. Mazumdar, Phys. Rev. D67, 046006 

(2003), hep-th/0210254. 
[77] C. Germani and A. Kehagias (2009), 0908.0001. 
[78] C. Armendariz-Picon, T. Damour, and V. F. Mukhanov, 

Phys. Lett. B458, 209 (1999), hep-th/9904075. 
[79] A. Gruzinov (2004), astro-ph/0401520. 
[80] M. S. Turner, Phys. Rev. D 28, 1243 (1983). 
[81] V. F. Mukhanov, H. A. Feldman, and R. H. Branden- 

berger, Phys. Rept. 215, 203 (1992). 
[82] H. Ruegg and M. Ruiz-Ahaba, Int. J. Mod. Phys. A19, 

3265 (2004), hep-th/0304245. 
[83] T. Koivisto and D. F. Mota, Phys. Rev. D75, 023518 

(2007), hep-th/0609155. 
[84] T. Koivisto, Phys. Rev. D77, 123513 (2008), 0803.3399. 
[85] J. B. Jimenez and A. L. Maroto, JCAP 0903, 016 (2009), 

0811.0566. 

[86] J. B. Jimenez, T. S. Koivisto, A. L. Maroto, and D. F. 

Mota (2009), 0907.3648. 
[87] M. Novello and J. M. Salim, Phys. Rev. D20, 377 (1979). 
[88] M. Novello, E. Goulart, J. M. Salim, and S. E. 

Perez Bergliaffa, Class. Quant. Grav. 24, 3021 (2007), 

gr-qc/0610043. 

[89] I. T. Drummond and S. J. Hathrell, Phys. Rev. D22, 343 
(1980). 

[90] K. Bamba and S. D. Odintsov, JCAP 0804, 024 (2008), 
0801.0954. 

[91] J. Sadeghi, M. R. Setare, and A. Banijamah (2009), 
0906.0713. 

[92] K. Bamba and S. Nojiri (2008), 0811.0150. 

[93] C. M. WiU and J. Nordtvedt, Kenneth, Astrophys. J. 

177, 757 (1972). 
[94] T. Eguchi, P. B. Gilkey, and A. J. Hanson, Phys. Rept. 

66, 213 (1980). 
[95] See figure 3.1 of Ref. [9^. 

[96] The nontrivial case of scalar potential including both 
quartic and quadratic terms has been solved just recently. 



